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Abstract. The conventional picture of disk accretion onto magnetized neutron stars has been challenged by the 
spin changes observed in a few X-ray pulsars, and by theoretical results from numerical simulations of disk- 
magnetized star interactions. These indicate possible accretion during the propeller regime and the spin-down 
torque increasing with the accretion rate. Here we present a model for the accretion torque exerted by the disk on 
a magnetized neutron star, assuming accretion continues even for rapid rotators. The accretion torque is shown 
to have some different characteristics from that in the conventional model, but in accord with observations and 
numerical calculations of accretion-powered magnetized neutron stars. We also discuss its possible applications 
to the spin evolution in X-ray pulsars. 
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1. Introduction 

The interaction between a magnetized, rotating star and 
a surrounding accretion disk takes place in a variety of as- 
trophysical systems, including T Tauri stars, cataclysmic 
variables, and neutron star X-ray binaries (Frank, King, 
& Raine 2002). Magnetic fields play an important role 
in transferring angular momentum between the accreting 
star and the disk. A detailed model for the disk-star inter- 
action was developed by Ghosh & Lamb (|1979al fl979b). In 
this model the stellar magnetic field lines are assumed to 
penetrate the accretion disk, and become twisted because 
of the differential rotation between the star and the disk. 
The resulting torque can spin up or down the central star, 
depending on the star's rotation, magnetic field strength, 
and mass accretion rate. The total torque N exerted on a 
star with mass M contains two components, 

N = N + AT mag . (1) 

The torque Nq carried by the material falling onto the star 
is generally taken to be 

A^o = MR 2 n (2) 

where Slo — {GM / R^) 1 / 2 is the (Keplerian) angular veloc- 
ity at the inner radius i?o of the disk, G the gravitational 
constant, and M the accretion rate. The magnetic torque 
-^mag results from the azimuthal field component gen- 
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erated by shear motion between the disk and the vertical 
field component -Bz 1 , 

/>oo 

N mag = / R 2 B 4> B z dR. (3) 

The torque N mllg can be positive and negative, depending 
on the value of the "fastness parameter" u> = f2*/fio = 
(Ra/R c ) 3 / 2 , where f2» is the angular velocity of the star 
and R c = (GM/Q 2 ) 1 / 3 the corotation radius. A rapidly 
rotating star (uj ^ 1) can be spun down with steady ac- 
cretion. The Ghosh & Lamb picture has been modified 
and extended in several directions by, for example, Konigl 
(EMI, Campbell ifT^jl . Lovelace et al. lP35|) . Wang 
(IUSI), Yi (P95|> . Li & Wang <HH|), Li, Wickramasinghe, 
& Ritdiger l(T$M}> . 

However, the above standard model has been chal- 
lenged by recent observations on the spin changes in sev- 
eral X-ray pulsars (Bildsten et al. 12197) and accreting mil- 
lisecond pulsars (Galloway et al. 2002), and by theoretical 
results from 2- and 3-dimensional numerical simulations of 
disk accretion to magnetized, rotating stars (Romanova et 
al. I2UU21 13UU51 12004[> . These observational and theoretical 
facts have raised a number of puzzling issues, including 
(1) abrupt spin reversals in X-ray pulsars, (2) spin-down 
in millisecond X-ray pulsars throughout the outburst, (3) 
spin-down rates increasing with accretion rate, and (4) ac- 
cretion in the "propeller" regime. These have motivated 

1 Here we have adopted cylindrical coordinate (R, <f>, z) cen- 
tered on the star, and the disk is assumed to be located on 
the z = plane, perpendicular to the star's spin and magnetic 
axes. 
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more recent investigations on the torques exerted by an 
accretion disk on a magnetized star (Nelson et al. 1997 
Yi, Wheeler, & Vishniac 1TW7I Torkelsson lTOOSI Lovelace, 
Romanova, & Bisnovatvi-Kogan 119991 Locsei & Melatos 
13001 Rappaport et al. 13004) . 

In this paper, attempting to solve some of the puz- 
zles related to disk-accreting neutron stars, we have con- 
structed a new model for the accretion torque based on 
some modifications of the Ghosh & Lamb picture. We in- 
troduce the model in section 2, and apply it to spin evolu- 
tion in X-ray pulsars in section 3. In section 4 we discuss 
the similarities and differences between this work and the 
competent model proposed by Lovelace et al. ©99 ). 



2. Model 

Here we adopt the typical assumptions. The central star 
contains a rotation-axis aligned dipolar magnetic field. At 
a cylindrical radius R from the star, the vertical compo- 
nent of the field B Z (R) = B*(R*/R) 3 7 where B* is the 
stellar surface field strength. This field is anchored in the 
stellar surface and also threads a thin-Keplerian accre- 
tion disk. The region above and below the disk contains 
low-density gas and the system exists in a steady-state. 
Assume that the magnetosphere is nearly force-free and 
reconnection takes place outside the disk (Wang ©95) , 
the azimuthal component of the field on the surface of 
the disk generated by rotation shear is given by 



motion of the disk material (Campbell 1992, Wang 1995), 



B,p 

~bI 



(n.-n K ) 
(n.-n K ) 



Cl* < Ok 
O* > Ok, 



(4) 



where the parameter 7 is of order unity, and Ok = Ok {R) 
is the Keplerian angular velocity in the disk. This form 
guarantees that the azimuthal pitch \B^/B Z \ cannot be 
larger than 7 over significant distances, for reasons of equi- 
librium and stability of the field above and below the disk 
plane fAlv 11351 113551 Uzdensky, Konigl, & Litwin 12002) ) . 

We assume that Eq. (4) applies in both accretion 
[u> < 1) and "soft" propeller (u> 1) cases. It is con- 
ventionally thought that, when M drops sufficiently, an 
accreting neutron star will transit from an accretor to a 
propeller, in which mass ejection occurs. However, this 
point of view is not fully consistent with recent observa- 
tions of transient, accretion-driven X-ray pulsars and with 
the energetics considerations (see Rappaport et al. 2004 
for a discussion). In stead, Spruit & Taam (1993) argued 
that even when uj £ 1 accretion cannot stop, and the disk 
may find a new accretion state with the inner disk ra- 
dius close to the corotation radius R c . Rappaport et al. 
(2004) further developed this idea and applied the results 
to compute the torques expected during the outbursts of 
the transient millisecond pulsars. 

The inner edge Rq of an magnetically truncated accre- 
tion disk is usually assumed to be determined by the con- 
dition that the magnetic stresses dominate over viscous 
stresses in the disk and begin to disrupt the Keplerian 



R -Af d(iZ2n) l 

.K o .«0O-DzO — M — |o, 



(5) 



where the subscript denotes quantities evaluated at R . 
The limitation of this definition is obvious if one combines 
Eqs. (4) and (5), and takes the limit as O approaches 
Oo: the left hand side of Eq. (5) becomes zero while the 
right hand side keeps to be a certain value, suggesting 
that it is not possible to calculate Rq when ui — > 1 (Li 
& Wickramasinghe 1997). However, the critical fastness 
parameter lo c beyond which spin-down occurs is always 
close to 1 (Wang Elil H2I2 Li & WanglJJiJ. Since in 
our model accretion occurs with u> extending from to 
1, we in stead take Rq to be at the magnetosphcric 
radius i? m determined by equating the ram pressure of 
the accreting matter with the magnetic pressure due to 
the dipole field of the neutron star (Davidson & Ostriker 
[TOTS! Lamb, Pethick, & Pines 11575^ . 



Ra 



Rm — ( 



I' 



2GMM 2 



1/7 



(6) 



where /1 = B*Rl is the magnetic moment of the neutron 
star. Observational evidence supporting this expression 
has been found in the transient X-ray pulsar A0535+26 
(Li 1997). We have not taken Rq = R c for uj > 1 as in 
Rappaport et al. (20Q3J). The reason is that, if the disk 
has already been disrupted at R m (> R c ), the accretion 
flow may pile up around the magnetosphere (Romanova 
et al. l2004f) . it is unknown whether the accretion flow can 
still keep the "disk" structure between R c and R m . 

We now estimate the torque exerted on the neutron 
star by the surrounding accretion disk. When the disk 
rotation is Keplerian, we integrate the magnetic torque in 
Eq. (3) over the disk from Rq to infinity combining with 
Eq. (4), and get the following expression, 



(^ )( 1-2 W + 2|!), W <1 



mag 



UJ > 1. 



(7) 



For the material torque Nq, Eq. (2) implies 100% effi- 
cient angular momentum transfer. It actually results from 
the magnetospheric torque arising from the shearing mo- 
tion between the corotating magnetosphere and the non- 
Keplerian boundary layer in the disk, which connects the 
magnetosphere and the outer Keplerian disk. The multi- 
dimensional calculations by Romanova et al. ( 2002 , 2004 ) 
indicate that most of the angular momentum carried by 
the accreting material in the magnetosphere has been 
transferred to the magnetic field. From the investigations 
of magnetized accretion disks (e.g., Ghosh & Lamb l*H?79al 
Campbell 1X3571 Erkut & Alpar l2"00ljl we find that the an- 
gular velocity in the boundary layer can be roughly rep- 
resented with the following form 



Q(R) = n (R/R y 



(8) 
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and the magnetospheric torque can be evaluated according 
to Eq. (4), 

N = - R 2 B^B z dR 
JRo-A 
rRo 



p -tto (~) p 

/ i?- 4 [i-(Jf)(#r i ]di? 

JR -A "0 -n-0 



(V2j8)M{GMR ) 1/2 {l - uj) 



(9) 



where A is the width of the boundary layer and 5 = 
A/R < 1 (Ghosh & Lamb HSBSi) . 

Obviously the real form of Q (R) must be more com- 
plicated than Eq. (8). However, Eq. (9) presents a reason- 
able order-of-magnitude estimate of the magnetospheric 
torque. Thus we let the torque Nq take the following form 
in both accretion and propeller cases (see also Menou et 
al. - 



(10) 



iV = £,M(GMR ) 1/2 (1 - uj), 



where the parameter £ = \/2^8 depends on the struc- 
ture of the magnetosphere and < £ ;$ 1. This implies 
that the accreting material in the boundary layer may 
lose some fraction of the angular momentum, probably, 
through the viscous torque and/or the magnetic centrifugal 
winds / outflows (e.g. Erkut & Alpar[2004). In the present 
analysis we always take £ = 1, and Eq. (10) recovers to 
Eq. (2) when The total torque exerted on the star 

by the disk is then 



N = N + N n 



M(GMi? ) 1/2 K(l-^) + 4 2 (l 



M(GAfi? ) 1/2 [£(l -u) 



(: 



2w + ^-)], uj 

- 1)], UJ > 1. 



Figures la and lb show the dependence of the dimen- 
sionless torque N/M(GMR ) 1 / 2 and N/(n 2 /V2R 3 c ) on uj 
in solid curves with £ = 7 = 1, respectively. The dashed 
curves in the figures represent the results in the disk model 
of Wang ( 1993J) for comparison. In Fig. la there is no sin- 
gularity problem in this work when uj — ► 1, the parameter 
space range for uj for spin-down with accretion is also much 
wider. The critical fastness parameter lo c , at which N — 0, 
is uj c ~ 0.884. Furthermore, it indicates that, for constant 
M, the spin-down torque increases with £1* (oc uj), in line 
with the numerical simulation results by Romanova et al. 

Figure lb shows how the torque varies with uj oc M -3 / 7 
when (i and are invariant. One can see that the spin- 
down torque is not a monotonous function of uj (or M) 
when uj > ui c . The spin-down torque takes the maximum 
value when uj ~ 1.634. Thus a given spin-down torque can 
correspond to two values of uj. When uj > 1.634, the spin- 
down torque increases with M, which is opposite to the 
dependence in the standard model, but consistent both 
with the observational fact in the X-ray pulsar GX 1+4 



that the X-ray flux appears to be increasing with the spin- 
down torque (Chakrabarty 1995 ), and with the numerical 
simulation results by Romanova et al. {2004). 

Finally we point out that the derived results in this 
simplified model should be regarded as qualitative, and 
one should not take the values above very seriously. 

3. Applications to X-ray pulsars 

The long-term, continuous monitoring of accreting X- 
ray pulsars by BATSE on board Compton Gamma- Ray 
Observatory has provided new insight into the spin evo- 
lution of these systems (Bildsten et al. I1997fl . For the X- 
ray pulsar Cen X— 3, the BATSE data show that it ex- 
hibits 10—100 day intervals of steady spin-up and spin- 
down, with a transition timescale 10 days. Similar bi- 
modal torque (or torque reversal) behavior has been ob- 
served in other pulsars, including 4U 1626—67, GX 1+4, 
and OAO 1657—415. These pulsars seem to be subject to 
instantaneous torques of roughly comparable magnitude 
iS M(GMRc) 1 / 2 , and differentiate themselves only by the 
timescale for reversals of sign. It is hard for the standard 
Ghosh & Lamb type model to explain a sudden torque 
reversal with nearly constant \N\ unless the M variation 
is discontinuous and fine tuned. 

Many mechanisms have been proposed to explain the 
torque reversals in accretion-powered X-ray pulsars. They 
can be roughly divided into two categories, one for the 
change of the direction of disk rotation, the other for bi- 
modal magnetic torque. In the former type of models, 
van Kerkwijk et al. (1998) suggested that the accretion 
disk may be subject to a warping instability because of 
the irradiation from the neutron star. The inner part of 
the accretion disk will flip over and rotate in the oppo- 
site direction, which would lead to a torque reversal. The 
spin-down can also be explained if the disk rotation is 
retrograde (Makishima et al. 119881 Nelson et al. I1997|l . 
In the latter, Torkelsson Q1998[) suggested that the mag- 
netic torque reversals may be explained by the presence 
of the intrinsic magnetic field in the accretion disk. The 
orientation of the disk field is not determined by the dif- 
ference between the angular velocities of the star and of 
the disk but is rather a free parameter. Thus the direction 
of the magnetic torque between the two is arbitrary. A bi- 
modal magnetic torque was also proposed by Lovelace et 
al. (1999, sec below). Locsei & Melatos (2004) presented a 
model of the disk-magnetosphere interaction, which adds 
diffusion of the stellar magnetic field to the disk. In cer- 
tain conditions, the system possesses two stable equilibria, 
corresponding to spin-up and spin-down. 

An alternative model for the transitions between spin- 
up and spin-down was suggested by Yi et al. (1997), 
in which the torque reversals are caused by alterna- 
tion between a Keplerian thin disk and a sub-Keplerian, 
advection-dominated accretion flow (ADAF) with small 
changes in the accretion rate. When M becomes smaller 
than a critical value M cr , the inner part of the accretion 
disk may make a transition from a primarily Keplerian 
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flow to a substantially sub-Keplerian ADAF, in which the 
angular velocity Sl'(R) = A&Ijz(R) with A < 1 (Narayan 
& Yi I1995[l . In this case the corotation radius becomes 
R' c = A 2 / 3 R C (Here we use the prime to denote quantities 
in ADAF). The dynamical changes in the disk structure 
lead to the slow (lu < u> c , spin-up) and rapid (u/ > lo' c , 
spin-down) rotator stage alternatively with M ~ M cr . 
This scenario is attractive because it is consistent with 
the transition time scales - the standard disk- ADAF tran- 
sition will occur on a thermal time t t h ~ (aSlo) -1 ~ 10 3 
s, while the interval between torque transition is set by 
the changes in M on a global viscous time scale. However, 
some issues in this model suggest that it needs to be im- 
proved. First, at equilibrium spin the critical fastness pa- 
rameter uj c — lj' c ~ 0.875 (see also Wang 1995). The nar- 
row range of w' (~ 0.875 — 1) for spin-down requires that 
the physical parameters (e.g., M, B*) should still be fine- 
tuned for various X-ray pulsars 2 . Second, as the Ghosh & 
Lamb model, this model also predicts decreased spin-down 
rate when M increases. 

In the present work, we have assumed that accretion 
can continue even when u) it 1- During the disk transition 
the inner radius of the disk may change little, i.e., R' Q ~ 
Ro, since the total (thermal and kinetic) energy density 
in the disk remains to be unchanged (Yi et al. 1997). The 
torque in the ADAF case can be derived in the similar 
way as in the Keplerian disk case, 



N' 



AM(GMR' Y/ 2 [^1 - u') + #(1 - 2uj' + %£■)], 
AM(GMR' )^iai «>') + #(g§r ~ 1)], J > 1 



where to' = (R^/R'J 3 / 2 = uj/A. Note that the critical fast- 
ness parameter uj' c (for N' = 0) equals uj c , but the equi- 
librium period would become longer by a factor of X/A, 
and the system begins to evolve towards the newly deter- 
mined equilibrium after the transition. Figure 2 shows the 
torque N and N' in units of M(GMR C ) X / 2 as a function 
of oj with £ = 7=1. For u between ~ 0.2 and ~ 0.8, if M 
is around M cr , small changes in M can cause the torque 
reversal behavior with comparable spin-up/down torques 
< M(GMRc) 1/2 . 

We have applied the model to explain the spin change 
in several X-ray pulsars. As an illustration, Fig. 3 shows 
the fit to the observed spin reversal in GX 1+4. Similar 
as in Yi et al. dH27|, we assume that M keeps a linear 
increase or decrease as an approximation to more complex 
M variations on longer timescales. The transition occurs 
on a timescale <C |f2*/f2*| before and after the transition. 
In Fig. 3, the observed torque reversal event is reproduced 
by taking B* = 1.0 x 10 13 G, M = 6.0 x 10 15 gs"\ 
dM/dt = 3.25 x 10 15 gs^ 1 yr" 1 , and A = 0.3 (we take 
the neutron star mass M = 1AM®, and radius R* = 10 6 
cm). We note, however, that our estimated values (of B* 
and M) can always be rescaled by changing A and 7. The 
critical accretion rate for the transition is assumed to be 
M cr = 6.5 x 10 15 gs- 1 . 

2 For oblique rotators, there could be no spin-down allowed 
with ui or u' < 1 (Wang ESSE). 



4. Discussion 

Assuming accretion continues during the (soft) propeller 
stage, we have constructed a model for the magnetized 
torque on accreting neutron stars, which changes continu- 
ously when co varies from to it 1 (without the singular- 
ity problem), and the spin-down torque increases with M. 
Following Yi et al. (199ZJ), we have applied the model to 
explain the torque reversal phenomena observed in a few 
X-ray pulsars. 

The same problem has been investigated by Lovelace 
et al. I|1999|l . These authors have developed a model for 
magnetic "propeller" -driven outflows for rapidly rotating 
neutron stars accreting from a disk. They find that the 
inner radius Rq of the disk depends not only on /1 and M, 
but also on the star's rotation rate 17*. The most interest- 
ing result is that for a given value of f2*, there may exist 
two solutions of Ro, with one > R c and the other < R c , 
corresponding to two possible equilibrium configurations. 
In a transition of Rq, the propeller varies between being 
"off' and being "on" , leading to jumps between spin-up 
and spin-down. The transitions are assumed to be stochas- 
tic or chaotic in nature, and could be triggered by small 
variations in the accretion flow and magnetic field config- 
uration. The ratio of spin-down to spin-up torque is also 
shown to be of order unity. 

This model shows some similar features as our work. 
For example, in both models (1) discontinuous transitions 
ra ^g.^ magnetic torque are invoked to account for the 
observed abrupt spin changes in X-ray pulsars, (2) the 
torques are of comparable magnitude, and (3) the spin- 
down torques increases with M. Especially, Eq. (10) sug- 
gests loss of angular momentum at the inner edge of the 
disk, probably by the magnetic-driven outflows proposed 
by Lovelace et al. I|1999(l . 

However, there exist significant differences between the 
two works. First, in the propeller regime, Lovelace et al. 
(199!jJ) let most of the accreting material be ejected from 
the star, while we still allow mass accretion to guarantees 
no significant variations in the X-ray luminosities during 
the transition of spin evolution, implying that to 3> 1 in 
Lovelace et al. I|1999f) . and lo it 1 in this work. Second, the 
torque reversals in our model are caused by the change 
of the dynamical structure of the disk, as suggested by 
Yi et al. (1997), at a certain, critical mass accretion 
rate. In Lovelace et al. I)1999|l the transition is stochastic, 
and could take place at any luminosities. Third, we have 
adopted the traditional Alfven radius as the inner radius 
Ro of the disk in both accretion and propeller phases. In 
Lovelace et al. (1999) Ro depends on f2* also, and its value 
increases with f2* in the propeller case. Obviously the hy- 
pothesis of £l*-depending Ro is more realistic, but how Ro 
changes with has not been well understood 3 . Future 



3 Romanova et al. J2(jy4} indeed observed increasing Ro with 
SI* in their calculations. But this change is caused by the fact 
that the initial dipolar field becomes non-dipolar when the star 
rotates very rapidly. 
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observations are expected to provide more stringent con- 
straints on the mechanisms suggested in both models. 
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Fig. 1. The solid lines show the dimensionless torque (a) N / M(GMR ) 1 / 2 and (b) N/()i 2 /3R%) as a function of the 
fastness parameter w. The dashed lines represent the standard model. 
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Fig. 2. The dimensionless torque 7V/M(GMi? c ) 1 / 2 (solid 
line) and N' / M{GMR C ) 1 / 2 (dashed line) as a function of 
the fastness parameter u. 
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Fig. 3. Torque reversal event in GX 1+4. The points are 
the observed data and the solid line corresponds to the 
fitting result. 



